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f~* > ; , Abstract We present the results of a spatial correlation analysis between uniden- 

tified gamma-ray sources in the 3EG catalog and candidates to super- 
nova remnants. This work extends a previous study made by Romero, 
Benaglia & Torres (1999). We also consider the gamma-ray emission 
^ ' variability and the spectral index for the sources with positional coinci- 

'. dence. 
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1. SNRS AS GAMMA-RAY EMITTERS 

^ ' Supernova remnants (SNRs) are usually considered as the main sources 

^ . of cosmic rays with energies below ~ 10 15 eV. Both electrons and pro- 

tons are believed to be accelerated by Fermi mechanism in the expanding 
shock front of these objects. The electrons produce synchrotron emis- 
sion detected at radio wavelengths whereas the interactions between rel- 
ativistic protons and ambient nuclei can produce neutral pions, which 
quickly decay yielding 7-ray emission at energies E > lOOMeV, in the 
EGRET range. The expected 7-ray flux at Earth was given by Drury et 
al. (1994), 

F(> lOOMeV, -4.4x40-^ {^)^\ d-D 

where E$n is the energy of the SN in ergs, 6 is the fraction of the total 
energy of the explosion converted into cosmic ray energy, and n and 
d have their usual meaning, number density and distance, respectively. 
In most cases, the expected flux at GeV energies is far too low to be 
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detected in the range of EGRET, but the presence of nearby clouds can 
produce a significant enhancement of the 7-ray emission. Such scenario 
has been recently studied by Combi et al (1998) in relation with the 
source 3EG J1659-6251. 

As showed by Aharonian et al. (1994), when a SNR hits a cloud a 
part of the proton population can be transported into it by convection, 
and illuminate the cloud by subsequent p-p interaction, via pion decay. 
The resulting flux is 

F(> lOOMeV) ~ q -< Md ° u \ (1.2) 

where q y is the 7-ray emissivity per H-atom in the cloud. This parameter 
can be related with its value in the vicinity of Earth by g 7 = kq & , with 
k > 1. 

Other mechanisms, like relativistic bremsstrahlung and inverse Comp- 
ton losses, associated with the leptonic component, can also play a role if 
the electron density and/or the photon fields are high enough, see for in- 
stance Pohl (1996). Then, there are many theoretical reasons to expect 
a positional correlation between SNRs and unidentified 7-ray sources. 
We devote the rest of this work to analyze this correlation. 

2. POSITIONAL CORRELATION 

Possible physical correlation between SNRs and unidentified EGRET 
sources, on the basis of two dimensional positional coincidence, has been 
proposed since the release of COS-B data (Montmerle 1979) and the first 
EGRET (1EG) catalog. Sturner & Dermer (1995) suggested that some 
of the unidentified sources lying at low galactic latitudes \b\ < 10° might 
be associated with SNRs: of 37 detections, 13 overlapped SNR positions 
in the 1EG catalog. However, the statistical significance was not too 
high as to provide a strong confidence. Using the 2EG catalog, Sturner 
et al. (1996) repeated the analysis, and showed that 95% confidence 
contours of 7 unidentified EGRET sources overlapped SNRs, some of 
them appearing to be in interaction with molecular clouds. Similar re- 
sults were independently reported by Esposito et al. (1996), considering 
only radio-bright SNRs, and Yadigaroglu and Romani (1997), although 
they did not assessed the overall chance probability of these 2EG-catalog 
findings. The evolution in the number of coincidences and SNRs con- 
sidered in the different studies is shown in Table 1. Note that from the 
First to the Second EGRET catalogs, fourteen unidentified sources were 
discarded. In the Third catalog, 6 unidentified sources are possibly arti- 
facts produced by the strong emission of the Vela pulsar: these sources 
disappear in a map where the pulsed Vela emission is suppressed. 1 
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In Table 2 we show the 3EG sources that are positionally coincident 
with SNRs listed in the latest version of Green's catalog. From left to 
right we provide the 7-ray source name, the measured flux, the photon 
spectral index T, the variability index I (see below), the SNR identifica- 
tion, the angular distance between the best 7-ray source position and the 
center of the remnant, the size of the remnant in arcminutes, the SNR 
type (S for shell, F for filled-centre, and C for composite), and other 
positional coincidences found in our previous study (Romero, Benaglia 
and Torres 1999). It is interesting to note that in the 3EG catalog, not 
all the positional coincidences with SNRs are SNOBs, as was the case re- 
ported by Montmerle (1979) and Yadigaroglu and Romani (1997) using 
previous samples. 

We also consider whether some of the sources in our sample may be 
associated with recently proposed candidates to supernova remnants, 
presently not catalogued by Green (1998). Our interest in this search 
resides in the fact that young stellar objects, like recently formed black 
holes and pulsars, can still be associated with the gaseous remnant of the 
original supernova that created them. The diffuse non-thermal emission 
of the galactic disk, originated in the interaction of the leptonic compo- 
nent of the cosmic rays with the galactic magnetic field, is surely veiling 
many remnants of low surface brightness. Recent observational studies 
using filtering techniques in the analysis of large-scale radio data have 
revealed several new SNR candidates that are not yet included in the 
latest issue of Green's catalog. In general, these new candidates are 
much more extended than those previously known. There are 101 of 
these weak non-thermal structures detected so far in the Galaxy. This 
number significantly extends Green's (1998) catalogue. The list of these 
new candidates, and the references from where they were compiled, can 
be obtained from the paper by Torres et al. (2000). 

We have found that only 7 gamma-ray sources in our sample are po- 
sitionally coincident with non-thermal radio structures. The positional 
coincidences thus obtained are shown in Table 3, where we provide sim- 
ilar information as that given in Table 2. 

In order to estimate the statistical significance of these coincidences, 
we have numerically simulated a large number of synthetic sets of EGRET 
sources using the code described in the paper by Romero, Benaglia 
and Torres (1999). The results of this study are shown in Table 4 
where we provide results for different samples of unidentified sources 
and SNRs catalogs: the original 81 unidentified EGRET sources, the 75 
high confidence ones, and the 40 sources without any positional corre- 
lation with known galactic gamma-ray emitters (as reported by Romero 
et al. (2000)). For the latter sample we are interested in checking 
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Table 1. Evolution in the number of positional coincidences between SNRs and unidentified EGRET sources 



Catalog 
EGRET 


Unidentified 
detections 


Real 
coincidences 


Number of SNRs 
in Green's catalog 


Significance 
(statistical) 


lEG a 


37 


13(35%) 


182 


1.8a 


2EG b 


32 


7(22%) 


194 


? 


2EG C 


32 


5(16%) 


14 g 


? 


2EG d 


33 


10(30%) 


194 


h 


3EG e 


81 


22(27%) 


220 


5.7a 


3EG f 


75 


19(25%) 


220 


4.8a 



Table 2. Positional coincidences with supernova remnants. 



7-source 


F 1 


r 


I 


SNR 


A6 


Size 


Type 


Other 


(3EG J) 


(10~ 8 ph cm~ 2 s" 1 ) 








(deg) 


(arcmin) 




coincidences 


0542+2610 


14.7+3.2 


2.67+0.22 


3.16 


G180.0-1.7 


2.04 


180 


S 




0617+2238 


51.4+3.5 


2.01+0.06 


1.68 


G189.1+3.0 


0.11 


45 


s 


OB 


0631+0642 


14.3+3.4 


2.06+0.15 


1.52 


G205.5+0.5 


1.97 


220 


s 




0634+0521 


15.0+3.5 


2.03+0.26 


1.02 


G205.5+0.5 


2.03 


220 


s 


Of/OB 


1013-5915 


33.4+6.0 


2.32+0.13 


1.63 


G284.3-1.8 


0.65 


24 


s 




1102-6103 


32.5+6.2 


2.47+0.21 


1.86 


G290.1-0.8 


0.12 


19 xl4 


s 


WR/OB 


1410-6147 


64.2+8.8 


2.12+0.14 


1.22 


G312.4-0.4 


0.23 




s 


OB/Of 


1639-4702 


53.2+8.7 


2.5+0.18 


1.95 


G337.8-0.1 


0.07 


9 x6 


s 


OB 










G338.1+0.4 


0.65 


15 


s 












G338.3+0.0 


0.57 


8 


s 




1714-3857 


43.6+6.5 


2.30+0.20 


2.17 


G348.5+0.0 


0.47 


10 


s 












G348.5+0.1 


0.50 


15 


s 




1734-3232 


40.3+6.7 




2.90 


G355.6+0.0 


0.16 


6 x8 


s 


OB 


1744-3011 


63.9+7.1 


2.17+0.08 


1.80 


G359.0-0.9 


0.41 


23 


s 












G359.1-0.5 


0.25 


24 


s 




1746-2851 


119.9+7.4 


1.70+0.07 


2.00 


GO.0+0.0 


0.12 


3.5 x2.5 


s 












GO.3+0.0 


0.19 


15.8 


s 




1800-2338 


61.3+6.7 


2.10+0.10 


1.60 


G6.4-0.1 


0.17 


42 


c 




1824-1514 


35.2+6.5 


2.19+0.18 


3.00 


G16.8-l.10 


0.43 


30 x24 




OB 


1837-0423 


19.1 


2.71+0.44 


5.41 


G27.8+0.6 


0.58 


50 x30 


F 




1856+0114 


67.5+8.6 


1.93+0.10 


2.92 


G34.7-0.4 


0.17 


35 x27 


s 




1903+0550 


62.1+8.9 


2.38+0.17 


2.28 


G39.2-0.3 


0.41 


8 x6 


s 




2016+3657 


34.7+5.7 


2.09+0.11 


2.06 


G74.9+1.2 


0.26 


8 x6 


F 


WR/OB 


2020+4017 


123.7+6.7 


2.08+0.04 


1.12 


G78.2+2.1 


0.15 


60 


s 


OB 



Table 3. Positional coincidences with candidates to SNRs. 



7-source 


F,< 


r 


I 


SNR 


A0 


Size 


Type 


(3EG J) 


(10~ 8 ph cm -2 s _1 ) 








(deg) 


(arcmin) 




0903-3531 


22.8+4.6 


2.66+0.24 


1.55 


260.20+1.40 


6.05 


1560 


S 


1631-4033 


12.7+3.0 


2.25+0.27 


1.26 


342.60+8.20 


3.11 


360 


S 


1638-5155 


67.3+14.2 


2.56+0.21 


2.44 


333.00+0.00 


3.50 


360 


s 


1704-4732 


117.8+27.1 


1.86+0.33 


2.96 


340.80-4.80 


1.22 


174 


s 


1717-2737 


17.8+3.6 


2.33+0.15 


2.00 


356.90+8.50 


2.66 


570 


s 


1834-2803 


16.2+2.9 


2.62+0.20 


2.83 


7.30-5.30 


3.91 


432 


s 


1837-0606 


49.6+7.9 


1.82+0.14 


2.41 


27.00+0.50 


1.14 


420 


s 



the possible association with candidates to SNRs. We conclude that 
there is strong statistical evidence suggesting that some 3EG unidenti- 
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Table 4. Statistical results. G stands for Green's SNRs and C for candidates. 



Number 
of SNRs 



Unidentified 
Sources 



coincidences coincidences probability 



Real Simulated Poisson 



321 (G + C) 
101 (C) 



220 (G) 
220 (G) 



81 
75 
75 
40 



22 7.8 ± 2.5 1.5 xl0~ 5 

19 7.0 ± 2.4 9.9 xl0~ 5 

30 22.4 ± 3.8 0.02 

7 10.4 ± 2.7 0.07 



fled EGRET sources must be physically associated with SNRs in Green's 
Catalog. At the same time, there is no statistical evidence suggesting 
that the 3EG sources analyzed in our sample are physically associated 
with candidates to SNRs: the number of real positional coincidences 
with them is totally compatible with, and even lower than, the result of 
a random association. This latter result appear to be the consequence 
of the candidates to SNRs being much more extended objects, thus im- 
proving the random coincidences with any given population, and does 
not discard in itself that some of the coincidences could be physical ones. 



3. SPECTRAL AND VARIABILITY INDICES 



In Figure 1 we present the distribution of the spectral index for all 
unidentified sources which resulted to have positional correlation with 
SNRs and candidates. All values are compatible with Fermi-like accel- 
eration processes that could happen in the strong shocks generated by 
the explosion of the supernova, or be the outcome of the interaction be- 
tween relativistic material and a nearby cloud, or -for the lowest ones- 
be associated with gamma-ray emission by pulsars. The distribution 
peaks around T = 2.1, but some sources have a large spectral index of 
about 2.8. Recall that the steepest measured spectral index for pulsars 
is around 2.2, so this is pointing against a pulsar origin for several of 
these gamma-ray sources. 

We now assess the possible long term variability of the sources. We 
define a mean weighted value for the EGRET flux as: 



N vp is the number of viewing periods for each gamma-ray source. F(i) 
is the observed flux in the i th -period, whereas e(i) is the corresponding 




-i 



(1.3) 
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Gamma-ray spectral index 



Figure 1 Distribution of spectral index for the sources with positional coincidences 
with SNRs. The clearer boxes stand for the sources coincident with candidates to 
SNRs. 



error in the observed flux. 2 We then define the fluctuation index n as: 
fi = 100 x o" sc j x (F)~ . In this expression, cr s d is the standard deviation 
of the flux measurements. In order to remove as far as possible any 
spurious variability introduced by the observing system, we computed 
the fluctuation index /i for the confirmed gamma-ray pulsars in the 3EG 
catalog. The identification by Kuiper et al. (2000) is not included 
because of the blazar contamination of the EGRET flux, observed by 
these authors. We adopt the physical criterion that pulsars are non- 
variable gamma-ray sources. Then, any non-null /i-value for pulsars 
is attributed to experimental uncertainty. We then define an averaged 
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statistical index of variability, /, as 

T /^source 




(1.4) 



< M >pul! 



Isars 



In terms of the averaged index I, the adopted variability criterion is 
then that variable sources will be those with / > 2.5, which is 3a away 
from the statistical variability of pulsars. 3 

In Figure 2 we show the distribution for the variability index I of those 
sources that are positionally coincident with SNRs. We see that most 
of the 19 sources positionally related with Green's SNRs have I < 2.5, 
being their mean value 2.17. Indeed, 12 sources out of 19 have I < 2. 
There are three sources with very high /-index: 3EG J1824-1514, J1837- 
0423 and 3EG J0542+2610, with I = 3.00, I = 5.41 and I = 3.16, 
respectively. These sources show spectral indices of 2.19, 2.71 and 2.67. 
The source 3EG J1824-1514 has been recently proposed by Paredes et al. 
(2000) as a faint microquasar detected through VLBI observations. This 
and other sources with I > 2.5, are very unlikely physically associated 
with the SNRs in the usual sense, i.e. being pulsars or nearby clouds in 
interaction with the swept up material. Nevertheless, compact objects 
like accreting black holes, or even isolated Kerr-Newmann black holes 
(Punsly et al. 2000) are interesting possibilities, they surely requires 
additional analysis. 
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1. Notes corresponding to Table 1 are as follow: a: Sturner and Dermer (1995). b: 
Sturner, Dermer and Mattox (1996). c: Esposito, Hunter, Kanbach and Sreekumar (1996). 
d: Yadigaroglu and Romani (1997). e: Romero, Benaglia and Torres (1999). f: This work, 
g: Only radio-bright SNRs, flux at 1 GHz greater than 100 Jy, were used, h: Computed for 
pairs. 



2. For those observations in which the significance (VTS in the EGRET catalog) is greater 
than 3<r, we took the error as e(i) = F(i)/\/TS. However, many of the observations are in 
fact upper bounds on the flux, with significance below 3<7. For these ones, we assume both 
F(i) and e(i) as half the value of the upper bound. 

3. Detailed comments on how this variability index compares with others, especially with 
the analysis made in W. Tompkins' Ph.D. thesis, are given elsewhere (Torres et al. 2000). 
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Figure 2 Distribution of the variability index for the sources with positional coin- 
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